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Abstract 

We present analytically the exact energy bound-states solutions of the 
Schrodinger equation in D-dimensions for an alternative (often used) pseudo- 
Coulomb potential-plus- ring-shaped potential of the form V(r) = — - + ^ + 
^> c °%^ + c by means of the conventional Nikiforov-Uvarov method. We give 

r sin J ° 

a clear recipe of how to obtain an explicit solution to the radial and angular 
parts of the wave functions in terms of orthogonal polynomials. The total 
energy of the system is different from the pseudo- Coulomb potential because 
of the contribution of the angular part. The general results obtained in this 
work can be reduced to the standard forms given in literature. 
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I. INTRODUCTION 



The important task of quantum mechanics is to find the exact bound-states solution of the 
Schrodinger equation for certain potentials of physical interest. Generally speaking, there 
are a few main traditional methods to study the exact solutions of quantum systems like 
the Coulomb, the harmonic oscillator [1,2], the pseudoharmonic [3,4] and the Kratzer [4,5] 
potentials. Additionally, in order to obtain the bound-states solutions of central potentials, 
one has to resort to numerical techniques or approximation schemes. For many of the 
quantum mechanical systems, most popular approximation methods such as shifted 1/N 
expansion [6], perturbation theory [7], path integral solution [8], algebraic methods with 
the SUSYquantum mechanics method and the idea of shape invariance, further closely with 
the factorization mrthod [9], exact quantization rule [10,11], the conventional Nikiforov 
and Uvarov (NU) method [12-25]. Some of these methods have drawbacks in applications. 
Although some other methods give simple relations for the eigenvalues, however, they lead 
to a very complicated relations for the eigenf unctions. 

The study of exact solutions of the Schrodinger equation for a class of non-central po- 
tentials with a vector potential and a non-central scalar potential is of considerable in- 
terest in quantum chemistry [26-35]. In recent years, numerous studies [36-40] have been 
made in analyzing the bound states of an electron in a Coulomb field with simultaneous 
presence of Aharanov-Bohm (AB) [41] field, and/or a magnetic Dirac monopole [42], and 
Aharanov-Bohm plus oscillator (ABO) systems. In most of these studies, the eigenvalues 
and eigenfunctions are obtained by means of seperation of variables in spherical or other 
orthogonal curvilinear coordinate systems. The path integral for particles moving in non- 
central potentials is evaluated to derive the energy spectrum of this system analytically [43]. 
In addition, the idea of SUSY and shape invariance is also used to obtain exact solutions of 
such non-central but seperable potentials [44]. Very recently, the conventional NU method 
has been used to give a clear recipe of how to obtain an explicit exact bound-states solutions 
for the energy eigenvalues and their corresponding wave functions in terms of orthogonal 
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polynomials for a class of non-central potentials [45]. 

Recently, Chen and Dong [46] found a new ring-shaped potential and obtained the exact 
solution of the Schrodinger equation for the Coulomb potential plus this new ring-shaped po- 
tential which has possible applications to ring-shaped organic molecules like cyclic polyenes 
and benzene. Recently, Cheng and Dai [47], proposed a new potential consisting from the 
modified Kratzer's potential [48] plus the new proposed ring-shaped potential in [46] . They 
have presented the energy eigenvalues for this proposed exactly-solvable non-central poten- 
tial in three dimensional (i.e., D = 3)-Schr6dinger equation through the NU method. The 
two quantum systems solved by Refs [46,47] are closely relevant to each other as they deal 
with a Coulombic field interaction except for a slight additional centrifugal potential barrier 
acts as a repulsive core which is for any arbitrary angular momentum £ prevents collapse of 
the system in any dimensional space due to the slight perturbation to the original angular 
momentum barrier. In a very recent work [49], we have given a clear recipe of how to ob- 
tain analytically the exact energy eigenvalues and the corresponding wave functions of the 
D-dimensional space Schrodinger equation with the proposed potential in [47] by means of 
the conventional NU method. 

The conventional Nikiforov-Uvarov (NU) method [12], which received much interest in 
recent years, has been introduced for solving Schrodinger equation [13-21], Klein-Gordon 
[22,23], Dirac [24] and Salpeter [25] equations. The purpose of this work is to follow parallel 
solution to [47] and give exact complete solutions to the D-dimensional Schrodinger equation 
with pseudo-Coulomb potential plus ring-shaped potential, a Coulombic-like potential with 
an additional centrifugal barrier, for any arbitrary ^'-states using the conventional Nikiforov- 
Uvarov method. For some appropriate choices of parameters, the obtained expressions of 
energy eigenvalues and their corresponding wave functions for this alternative form of the 
potential coincide with the D-dimensional results given in [49] and with the three-dimensions 
results given in [47] in the limiting case of D = 3. 

This work is organized as follows: in section II, we shall briefly introduce the basic 
concepts of the NU method. Section III is mainly devoted to the exact solution of the 
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Schrodinger equation in D-dimensions for this quantum system by means of the NU method. 
Finally, the relevant results are discussed in section IV. 



II. BASIC CONCEPTS OF THE METHOD 

The NU method is based on reducing the second-order differential equation to a gener- 
alized equation of hypergeometric type [12]. In this sense, the Schrodinger equation, after 
employing an appropriate coordinate transformation s = s(r), transforms to the following 
form: 

+ ^r/n(s) + ^MS) = 0, (1) 

where a(s) and <j(s) are polynomials, at most of second- degree, and f(s) is a first-degree 
polynomial. Using a wave function, ip n (s), of the simple ansatz: 

Vv»(s) = (f> n (s)y n (s), (2) 

reduces (1) into an equation of a hypergeometric type 

a(s)y^s) + r(s)y' n (s) + Xy n (s) = 0, (3) 

where 

r( s ) = t( s ) + 2tt(s), t'(s) < 0, (5) 

and A is a parameter defined as 

A = A n = -nr'(s) - " ^ ~ l) <r"{s), n = 0, 1, 2, .... (6) 

The polynomial r(s) with the parameter s and prime factors show the differentials at first 
degree be negative. It is worthwhile to note that A or A n are obtained from a particular 
solution of the form y(s) = y n (s) which is a polynomial of degree n. Further, the other 



part y n (s) of the wave function (2) is the hypergeometric-type function whose polynomial 
solutions are given by Rodrigues relation 

where B n is the normalization constant and the weight function p{s) must satisfy the con- 
dition [12] 

± w ( s ) = If\w(s), w(s) = a(s)p(s). (8) 
as a(s) 

The function n and the parameter A are defined as 



a'(s) - r(s) 

7r(s) = ± 



i 



\ = k + n'(s). (10) 

In principle, since ir(s) has to be a polynomial of degree at most one, the expression under 
the square root sign in (9) can be arranged to be the square of a polynomial of first degree 
[12]. This is possible only if its discriminant is zero. In this case, an equation for k is 
obtained. After solving this equation, the obtained values of k are substituted in (9). In 
addition, by comparing equations (6) and (10), we obtain the energy eigenvalues. 

III. EXACT SOLUTIONS OF THE QUANTUM SYSTEM WITH THE NU 

METHOD 

A. Seperating variables of the Schrodinger equation 

The pseudo- Coulomb potential plus ring-shaped potential in spherical coordinates is 
defined as [47] 

I7 , M a b n cos 2 6 

V(r, 9) = — + - + P- r -^ T - ft + c, 11 
r r z r z sin v 

where a, b, c and (3 are positive real constants. The potential in (11) coincides with the 
potential introduced by Cheng-Dai [47] if the parameters take their proper values and also 



turns to be the modified Kratzer potential in the limiting case of f3 — [48]. In fact the 
energy spectrum for this potential can be obtained directly by considering it as special case 
of the general non-central seperable potentials [45] . 

Our aim is to derive analytically the energy spectrum for a moving particle in the presence 
of a potential (11) in a very simple way. The .D-dimensional space Schrodinger equation in 
spherical polar coordinates written for potential (11) takes the form [1,6] 



2/i 



„D-1 



Or 



Or I r 2 \ sin 9 09 



09 



sin 2 9 Oif 2 h 2 sin 2 9 



il>(r,9,ip) 



+ 



a b ^ 
-- + ^ + c-E 



(12) 



where \x 



wave function separated as follows 



being the reduced mass of the two particles and ip(r, 9, ip) being the total 



ij nem (r, 9, tp) = R(r)Y e m (9, tp), R{r) = r-^/ 2 g(r), Y e m (9, tp) = H(0)*(<p). (13) 



On substituting equation (13) into (12) leads to a set of second-order differential equations: 

d 2 $(v?) 



dip 2 



+ m 2 $(y?) = 0, 



(14) 



1 did 
Md9 V dl 



m 



2/1/3 cos 2 9 
sin 2 9 ~ h 2 sin 2 9 



+ £(£ + D - 2) 



H{9) = 0, 



(15) 



1 d ( rD _ x d\ £(£ + D - 2) 



,D-1 



dr 



dr 



R(r) + | 



a b 
E+----C 



R(r) = 0. (16) 



The solution in (16) is periodic and must satisfy the period boundary condition &(p + 2ir) 
$(</?) from which we obtain 

1 



$ m (y?) 



exp(±imip), m — 0,1,2, 



(17) 



Further, equation (15) representing the angular wave equation takes the simple form 



d 2 H{9) cos9dH{9) 



d9 2 



+ 



sin 9 d9 



+ 



£(£ + D - 2) 



m 2 + (2fip/h 2 ) cos 2 9 
shi 2 ! 



H{9) = 0, 



(18) 



which has already been solved before in our previous work [49]. We present the final results 
for the angular part of the wave function as 



H m >{9) 



{21' + -m')\ 
\ 2(f + m')\ 



sm(6) m 'Pi m '> m '\cos9), 



(19) 



where 



n 



= -<I±M + i^'+l)2 + 4£'(D-3), 



(20) 



with 



f= JD-Q 1 



+ -\I(D- 2) 2 + 4(n + \/m 2 + 2,i/3/R')(n+l + v 'm 2 + 2/43/R ), (21) 



and 



m' 2 = m 2 + 2///3/r. 



(22) 



On the other hand, we are left to solve equation (16). After lengthy, but straightforward, 
calculations, Equation (16), representing the radial wave equation can be rewritten as [6]: 



d 2 g(r) 
dr 2 



+ 



2/i 



(E-c) + 



2/m 1 u+ (2fib/ti 



2\H 



h 2y 1 h 2 r 



g{r) = 0, 



(23) 



where 



v = _(M - 1)(M - 3), M = D + 21. 



(24) 



The two particles in equation (18) interacting via Coulombic-like field have a slight change 
in the angular momentum barrier acts as a repulsive core which for any arbitrary t prevents 
collapse of the system in any space dimension due to the additional centrifugal potential 
barrier. Thus, equation (23) has to be solved by using the NU method in the next subsection. 

B. The solutions of the radial equation 



The aim of this subsection is to solve the problem with a different radial separation 
function g(r) in any arbitrary dimensions. We now study the bound-states (real) solution 
E < c of equation (23). Letting 



2/i 

h 2 ' 



2 fia 



2[ib 



e = \I- T 2(E - c), a = 7 = i/ + 2 , 



(25) 



and substituting these expressions in equation (23), one gets 

w> + w^-a g(r) = 



(26) 



To apply the conventional NU-method, equation (26) is compared with (1) and the following 
expressions are obtained 



r(r) = 0, a(r) = r, a(r) = — e r + ar — 7. 
Substituting the above expressions into equation (9) gives 

7r(r) = - ± -\J^e 2 r 2 + 4(k - a)r + 47 + 1. 



(27) 



(28) 



According to this conventional method, the expression in the square root be the square of a 
polynomial. Thus, the two roots k can be readily obtained as 



k — a ± £^47 + 1. 



In view of that, we arrive at the following four possible functions of ir(r) : 



(29) 



7r(r) = < 



1 + 



(30) 



er + \y/¥i + lj for k x = a + ey/¥f+T, 

- er + iV 4 7 + l] for fci = a + ey/¥f+T, 
+ [er - |V 4 7 + l] for ^2 = a - ey/¥f+T, 

— er — + 1 for k 2 = a — e^A^y + 1. 

The correct value of ir(r) is chosen such that the function r(r) given by equation (5) will 
have negative derivative [12]. So we can select the physical values to be 



k = a — e^A^y + 1 and n(r) — - — er — - ^47 + 1 



(31) 



which yield 



r(r) = -2er + (1 + ^7+ 1). 



(32) 



Using equations (6) and (10), the following expressions for A are obtained, respectively, 

A = \ n = 2Ne, N = 0, 1, 2, (33) 



A = a-e(l + ^4 7 + l). (34) 
So we can obtain the energy eigenvalues as 

2ua 2 /h 2 , . 

E N = c- E . ~ (35) 



(2N + 1 + yJ(M - 1)(M - 3) + 8jjb/h 2 + 1 



where 

(M - 1)(M - 3) = Av = (D - 2) 2 + 4£'(£' + D - 2) - - 1, (36) 

with £' defined in (21). Therefore, the final energy spectra in equation (35) take the following 
Coulombic-like form [7] 

E N , = c- ^ a2 l % N > = o i 2 , ... (37) 
2 (N 1 ) 



where 

2N + \J (D - 2) 2 + W(£' + D - 2) + 8^(6 - [3)/k 2 + lj , (38) 



iV' = iV + L + l = - 



is simply obtained by means of substituting equation (36) into (35). 

(i) If a = 2D e r e , b = D e r 2 , c = D e and D = 3 (modified Kratzer potential-plus-ring- 
shaped potential), the equation (35), with the help of equation (36), is transformed into the 
following form 

p . n ^pyjn 2 

j 2N + 1 + J (2n + l) 2 + 4m 2 + 4(2n + l)^/m 2 + 2fi(3/h 2 + 8fiD e r 2 /h 2 \ 

(39) 

and it is consistent with Ref. [47]. 

(ii) If a = 2D e r e , b = D e r 2 , c = D e , D = 3 and (5 — (modified Kratzer potential), the 
equation (35) is transformed into the form 



K = D. (40) 

+ 2n + yl + U{t + 1) + 8/i&/ft 2 ) 

and it is consistent with Ref. [48]. 

(iii) If a = 2D e r e , b = D e r 2 e and c = D e , the equation (35) is transformed into the form 

&Nnm — L> e 9 , (41 ) 

[2N + 1 + \j (D - 2) 2 + 4f {£' + D-2) + 8/i(D e r 2 - (3)/h 2 ) 

and it is consistent with Ref. [49]. 

(iv) If a = Ze 2 , b = c = and D = 3, the equation (35) is transformed into the form 

2nZ 2 e 4 /h 2 



E 



Nnm 



27V + 1 + Jl + 4 



n + ^jm 2 + 2fi(3 /h 2 ^j (n + ym 2 + 2^/Ti 2 + l) - 2^jU 2 

(42) 



where m, n, N — 0, 1, 2, ... and it is consistent with Ref. [46]. 

Let us now turn attention to find the radial wavefunctions for this potential. Using r(r), 
7r(r) and a(r) in equations (4) and (8), we find 

0(r) = r (^+i+i)/2 e - £r , (43) 

p (r) = r V37+T e -2er_ (44) 

Then from equation (7), we obtain 



y n (r) = B n r-^e 2 -^ (^W*Fr e -*r) > (45) 

and the wave function g(r) can be written in the form of the generalized Laguerre polyno- 
mials as 

g(r) = C N , L r L+1 e- £r L 2 N L+1 (2er), (46) 

where 
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2 



(D - 2) 2 + W (£' + D-2) + 8/i(6 - (3)/h 2 - 1 



(47) 



Finally, the radial wave functions of the Schrodinger equation are obtained 



R(r) = C N , L r L ^ D '^e-^L 2 N L+1 (2er), 



(48) 



where 



h 2 N' 



(49) 



with N' is given in equation (38) and Cn,l is the normalization constant to be determined 

oo 

below. Using the normalization condition, / R 2 ^)^' 1 ^ = 1, and the orthogonality re- 

o 

lation of the generalized Laguerre polynomials, / z n+1 e~ z [L^(z)] 2 dz = ( 2 "+ ? ?+ 1 K"-+ ? ?) ! ^ we 
have 



(2e) 2L + 3 N\ 



(50) 



/N ' L \\ 2(N + L + 1)(N + 2L + 1)!' 

Letting D = 3 and a = a, the equation (48) transforms into equation (20) in Ref.[46] 
Therefore, we may express the normalized total wave functions as 

ip(r,6,<p) 



(2e) 2i+3 (2f + !)(£' -m')\N\ 



\ 4vr(f + m')\(N + L + 1)(N + 2L + 1)! 



L -( D -^ 2 exp(-er)L 2 N L+l (2er) 



x sin(^) m 'Pi m '' m,) (cos^) exp(±im<p). 



(51) 



IV. RESULTS AND CONCLUSIONS 

In this paper, the Schrodinger equation in any arbitrary dimensions has been solved for 
its exact bound-states with a a pseudo- Coulomb potential plus ring-shaped potential by 
means of a simple conventional NU method. The analytical expressions for the total energy 
levels of this system is found to be different from the results obtained for the modified 
Kratzer's potential in [48] and also more general than the one obtained recently in D- 
dimensions [49]. Hence, the noncentral potentials treated in [45] can be introduced as 

11 



perturbation to the central potentials by adjusting the strength of the coupling constant (3 
in the non-central part in terms of the coupling constants of the central part. In addition, 
the angular part, the radial part and then the total wave functions are also found. Thus, the 
Schrodinger equation with a new non-central but seperable potential has also been studied 
(cf. [45,49] and the references therein). This method is very simple and useful in solving 
other complicated systems analytically without given a restiction conditions on the solution 
of some quantum systems as the case in the other models. Finally, we point out that these 
exact results obtained for this new proposed form of the potential (11) may have some 
interesting applications in the study of different quantum mechanical systems, atomic and 
molecular physics. 

ACKNOWLEDGMENTS 

This research was partially supported by the Scientific and Technological Research Coun- 
cil of Turkey. S.M. Ikhdair wishes to dedicate this work to his family for their love and 
assistance. 



12 



REFERENCES 



[1] L.I. Schiff, Quantum Mechanics 3rd edn (McGraw-Hill Book Co., New York, 1968). 

[2] L.D. Landau and E.M. Lifshitz, Quantum Mechanics (Non-Relativistic Theory) 3rd edn 
(Pergamon, New York, 1977). 

[3] S.M. Ikhdair and R. Sever, J. Mol. Struct ure:THEOCHEM 806 (2007) 155. 

[4] S.M. Ikhdair and R. Sever, preprint quant-ph/0702052. 

[5] S.M. Ikhdair and R. Sever, to appear in the Int. J. Mod. Phys. E (preprint quant- 
ph/0611065). 

[6] S.M. Ikhdair and R. Sever, Z. Phys. C 56 (1992) 155; C 58 (1993) 153; D 28 (1993) 1; 
Hadronic J. 15 (1992) 389; Int. J. Mod. Phys. A 18 (2003) 4215; A 19 (2004) 1771; A 
20 (2005) 4035; A 20 (2005) 6509; A 21 (2006) 2191; A 21 (2006) 3989; A 21 (2006) 
6699; Int. J. Mod. Phys. E (in press) (preprint hep-ph/0504176); S. Ikhdair et al, Tr. J. 
Phys. 16 (1992) 510; 17 (1993) 474. 

[7] S.M. Ikhdair and R. Sever, Int. J. Mod. Phys. A 21 (2006) 6465; J. Math. Chem. (DOI 
10.1007/sl0910-006-9080-2); J. Math. Chem. (preprint quant-ph/0604078); J. Mol. 
Structure:THEOCHEM (DOI:10.1016/j.theochem.2007.01.019). 

[8] §. Erkoc, R. Sever, Phys. Rev. D 30 (1984) 2117. 

[9] F. Cooper, A. Khare and U. Sukhatme, Phys. Rep. 251 (1995) 267. 

[10] Z.Q. Ma and B.W. Xu, Europhys. Lett. 69 (2005) 685. 

[11] S.H. Dong, D. Morales and J. Garcia-Ravelo, Int. J. Mod. Phys. E 16 (2007) 189. 

[12] A.F. Nikiforov and V. B. Uvarov, Special Functions of Mathematical Physics 
(Birkhauser, Bassel, 1988). 

[13] S.M. Ikhdair and R. Sever, Int. J. Theor. Phys. (DOI 10.1007/sl0773-006-9317-7; J. 



13 



Math. Chem. (DOI 10.1007/sl0910-006-9115-8). 

[14] C. Berkdemir, A. Berkdemir and R. Sever, Phys. Rev. C 72 (2005) 027001; J. Phys. A: 
Math. Gen. 39 (2006) 13455; A. Berkdemir, C. Berkdemir and R. Sever, Mod. Phys. 
Lett. A 21 (2006) 2087. 

[15] H. Egrifes, D. Demirhan and F. Buyukkihc, Phys. Scr. 59 (1999) 90; 60 (1999) 195. 

[16] H. Egrifes, D. Demirhan and F. Buyiikkrhg, Theor. Chem. Acc. 98 (1997) 192. 

[17] 6. Ye§ilta§, M. §im§ek, R. Sever and C. Tezcan, Phys. Scr. 67 (2003) 472. 

[18] M. Znojil, Phys. Lett. A 264 (1999) 108. 

[19] M. Akta§ and R. Sever, Mod. Phys. Lett. A 19 (2004) 2871. 

[20] 6. Ye§ilta§, Phys. Scr. 75 (2007) 41. 

[21] F. Yasuk, A. Durmus. and I. Boztosun, J. Math. Phys. 47 (2006) 082302. 

[22] S.M. Ikhdair and R. Sever, Ann. Phys. (Leipzig) 16 [3] (2007) 218. 

[23] M. §imsek and H. Egrifes, J. Phys. A: Math. Gen. 37 (2004) 4379. 

[24] H. Egrifes and R. Sever, Phys. Lett. A 344 (2005)117. 

[25] S.M. Ikhdair and R. Sever, Int. J. Mod. Phys. E (in press) (preprint quant-ph/0605045) 

[26] M. Kibler and T. Negadi, Int. J. Quantum Chem. 26 (1984) 405. 

[27] i. SSkmen, Phys. Lett. 118A (1986) 249. 

[28] M. Kibler and P. Winternitz, J. Phys. A 20 (1987) 4097. 

[29] L.V. Lutsenko et at, Teor. Mat. Fiz. 83 (1990) 419. 

[30] H. Hartmann et al, Theor. Chim. Acta 24 (1972) 201. 

[31] H. Hartmann and D. Schuch, Int. J. Quantum Chem. 18 (1980) 125. 



14 



[32] M.V. Carpido-Bernido and C. C. Bernido, Phys. Lett. 134A (1989) 315. 
[33] M.V. Carpido-Bernido, J. Phys. A 24 (1991) 3013. 

[34] O. F. Gal'bert, Y. L. Granovskii and A. S. Zhedabov, Phys. Lett. A 153 (1991) 177. 

[35] C. Quesne, J. Phys. A 21 (1988) 3093. 

[36] M. Kibler and T. Negadi, Phys. Lett. A 124 (1987) 42. 

[37] A. Guha and S. Mukherjee, J. Math. Phys. 28 (1989) 840. 

[38] G. E. Draganescu, C. Campiogotto and M. Kibler, Phys. Lett. A 170 (1992) 339. 

[39] M. Kibler and C. Campiogotto, Phys. Lett. A 181 (1993) 1. 

[40] V. M. Villalba, Phys. Lett. A 193 (1994) 218. 

[41] Y. Aharonov and D. Bohm, Phys. Rev. 115 (1959) 485. 

[42] P. A. M. Dirac, Proc. R. Soc. London Ser. A 133 (1931) 60. 

[43] B.P Mandal, Int. J. Mod. Phys. A 15 (2000) 1225. 

[44] B. GSniil and i. Zorba, Phys. Lett. A 269 (2000) 83. 

[45] S.M. Ikhdair and R. Sever, to appear in the Int. J. Theor. Phys. (preprint quant- 
ph//0702186) 

[46] C.Y. Chen and S.H. Dong, Phys. Lett. A 335 (2005) 374. 
[47] Y.F. Cheng and T.Q. Dai, Phys. Scr. 75 (2007) 274. 

[48] C. Berkdemir, A. Berkdemir and J.G. Han, Chem. Phys. Lett. 417 (2006) 326. 
[49] S.M. Ikhdair and R. Sever, preprint quant-ph/0703008. 



15 



